Strong correlation effects and optical conductivity in electron doped cuprates 
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We demonstrate that most features ascribed to strong correlation effects in various spectroscopies 
of the cuprates are captured by a calculation of the self-energy incorporating effects of spin and 
charge fluctuations. The self energy is calculated over the full doping range of electron-doped 
cuprates from half filling to the overdoped system. The spectral function reveals four subbands, two 
widely split incoherent bands representing the remnant of the split Hubbard bands, and two addi- 
tional coherent, spin- and charge-dressed in-gap bands split by a spin-density wave, which collapses 
in the overdoped regime. The incoherent features persist to high doping, producing a remnant Mott 
gap in the optical spectra, while transitions between the in-gap states lead to pseudogap features in 
the mid-infrared. 
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PACS numbers: 74.25.Gz, 74.72.-h, 74.20.Mn,74.25.Jb 

A key issue in cuprate physics is to understand the 
routes through which cuprates evolve from a 'Mott' in- 
sulator at half filhng to a high- Tc superconductor at 
optimal doping. Underdoped cuprates show the pres- 
ence of a 'pseudogap' or a region of depleted density 
of states extending to several hundred meV's around 
the Fermi energy and thus are very unconventional, 
but the overdoped system behaves closer to a normal 
Fermi liquid and is thus amenable to a largely conven- 
tional description. Along these lines, evidence of small 
Fermi surface pockets [IJ has been reported in the un- 
derdoped regime, while in overdoped Tl2Ba2Cu06+5, a 
large, three-dimensional Fermi surface consistent with 
LDA band structure calculations [2| is observed. In sharp 
contrast, optical studies point to a far more complex and 
puzzling picture in that even in the overdoped case, an 
absorption peak characteristic of the Mott gap contin- 
ues to persist [si IJ] in the spectra, suggesting that dop- 
ing introduces new in-gap states in which the pseudogap 
physics resides, but that otherwise the Mott gap persists 
at all dopings. Here we show that a relatively simple and 
transparent model of the electronic self energy, where 
the quasiparticles are dressed with spin wave excitations, 
captures the key experimentally observed features of the 
remarkable doping evolution in optical and other spec- 
troscopies, including in particular the persistence of the 
Mott gap in the overdoped regime. Our self-energy is also 
in good accord with quantum Monte Carlo computations 
in the parameter regimes where the latter are available. 

We evaluate the self-energy E as a convolution over the 
green function G and the interaction W ~ U^X (includ- 
ing both spin and charge contributions) as [5|, 
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G'(fc-|-(f, a',a;„ -I- Wp)r(fc,g,u;„,Wp)Im[x'^'^ (q.LOp)], (1) 

where a is the spin index and the prime over the q 
summation means that the summation is restricted to 



the magnetic Brillouin zone. In the underdoped re- 
gion, the pseudogap is modeled by an antiferromagnetic 
(AFM) order parameter, resulting in G, x ^-^d E be- 
coming 2x2 tensorsJGJ. We define a total self-energy as 
E* = USti -f S, where fi is the Pauli matrix along the 
x— direction and US is the AFM gap defined below. The 
self-energy E* contains essentially two energy scales: (i) it 
gives rise to the SDW with an additional renormalization 
of the overall quasiparticle dispersions in the low energy 
region, and (ii) at higher energies it produces the Hub- 
bard bands. We use a modified self-consistent scheme, 
referred to as quasiparticle— GW^ (QP— GPF)-scheme in 
which G and W are calculated from an approximate self- 
energy Eo(cl') = USti -I- (l — Z~^) Lol, where the renor- 
malization factor Z is adjusted self-consistently to match 
the self-energy E* at low energy. [7i, 'S| 

We find that near optimal doping spin waves|7ll8l.l9Llld. 
ll| dress the quasiparticles into a coherent in-gap state, 
while the incoherent high-energy features are remnants of 
the upper and lower Hubbard bands (U/LHBs)[12]. With 
underdoping the in-gap state develops into a spin density 
wave (SDW) state which opens a gap between the upper 
and lower magnetic bands (U/LMBs). The model also 
describes the high-energy kink or the waterfall effect seen 
in the electronic dispersion'13] as the crossover between 
coherent and incoherent features. 

The present calculations are restricted to the electron- 
doped cuprates in order to avoid possible complica- 
tions of nanoscale phase separation. The self-energy 
Eg splits the LDA-band, <^j: (modelled by tight-binding 
(TB) parameters [7|) into renormalized UMB {v = +) 
andLMB (v = -): 
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where ^- == (^^r ± ^j:^^)/2. The AFM magnetization S 

ai Q — (tt, tt) is calculated self-consistently at each dop- 
ing, assuming a doping dependent [/.[1J| In the present 
formalism U is renormalized by Z. The doping depen- 
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FIG. 1: (color online) Calculated optical conductivity (solid 
lines) compared with experiment (symbols of same color) for 
dopings X — 0.0 to a; = 0.15 as taken from Ref. 3 and for 
X = 0.17 from Ref. iJ. Experimental data for x = 0.17 is 
subtracted by a background contribution to match the other 
data set. The dashed line gives the background contribution 
added to the theoretical spectrum at a; = 0. inset: Inferred 
doping dependence of the scattering rate. 



dency of U is chosen such that ZU reproduces the pseu- 
dogap in both angle-resolved photoemission spectroscopy 
(ARPES) and optical spectra, while x, S, and Z are de- 
termined self-consistently. 15| Finally, the vertex correc- 
tion r(fc, q, oj, ujp) in Eq. [T] is taken as its first order ap- 
proximation (Ward's identity) as T {k , q, ui , ujp) — l/Z. 
Since the fc-dependence of S is weak, we further simplify 
the calculation by assuming a fc- independent S, which 
we calculate at a representative point k = (7r/2, 7r/2). 

We begin by discussing the optical spectra of Fig. 1. 
The frequency dependent optical conductivity, (t(w), is 
calculated using standard linear response theory in the 
AFM state, [16[ for the full doping range from the half- 
filled state to the quantum critical point (QCP) in the 
overdoped region. A very good level of accord is seen 
with the experimental resultsjSJ, |j|. To fit the Drude 
conductivity, we have introduced an impurity scattering 
rate r which is found to have a strong doping dependence 
(inset). At high energy we include a doping-dependent 
background contribution, presumably associated with in- 
terband transitions to higher-lying bands not included in 
the present calculations. We have used the same energy- 
dependence of the background for all dopings with an 
intensity that decreases smoothly with doping. The red 
dashed line in Fig. 1 shows this interband contribution 
for the X = 0.0 spectrum. 

Interestingly, the spectra show a nearly isosbetic (equal 
absorption) point near 1.4 eV, consistent with the experi- 
mental behavior. The doping evolution is completely dif- 
ferent on opposite sides of this isosbetic point. Above this 



point the spectrum is dominated by a broad hump fea- 
ture above 1.5 eV, a signature of the Mott gap. At half- 
filling, only this feature is present and the optical spec- 
trum shows an insulating gap whose energy, structure, 
and intensity match remarkably with measurements [3|. 
As doping increases the high energy peak shifts to higher 
energy and broadens. 

Below the isosbectic point there is little spectral weight 
at half filling, but as doping increases spectral weight is 
gradually transferred from the higher energy region to 
the mid-infrared (MIR) one. The lower energy spectrum 
is associated with a Drude peak related to intraband 
transitions and a mid-infrared peak associated with tran- 
sitions across the magnetic gap [the pseudogap for the 
electron-doped cuprates]. With doping, this peak shifts 
to lower energy as the magnetic gap collapses and grad- 
ually sharpens due to the doping dependent scattering 
rate, see inset to Fig. 1. Note that at a; = 0.17, when 
the pseudogap has collapsed, Mott-gap features still per- 
sist in the spectrum. The present mean-field calculation 
overestimates the Neel temperature T^, but this can be 
corrected by including critical fluctuations, grave a la the 
Mermin- Wagner theorem [17.]. 

The origin of these features can be understood by look- 
ing at the doping dependence of the momentum-resolved 
spectral weight in Fig. 2. In the overdoped case in 
Fig. 2(e), a kink due to the bosonic coupling reproduces 
the waterfall effect below Ep , with a corresponding effect 
above Ep, splitting the spectrum into an effective three- 
band behavior, with UHB, LHB, and in-gap states. The 
features in the optical spectra are associated with transi- 
tions between these bands: The residual Mott gap arises 
from the transition from the LMB to the incoherent UHB 
[or from LHB to UMB], while the Drude term is asso- 
ciated with intraband transitions near the Fermi level. 
At lower doping an AFM gap opens in the coherent in- 
gap statesflS'l, leading to the UMB/LMB splitting and a 
four band behavior similar to that seen in QMC cluster 
calculations [l9J- Consistent with the QMC calculations, 
the coherent in-gap bands are dressed by magnetic quasi- 
particles. As the magnetic gap opens, the MIR feature 
in the optical spectra, being associated with transitions 
across this gap, shifts to higher energy. 

The quality of self-consistency of our scheme can be 
assessed by noting that the final coherent bands have 
nearly the same dispersion as the Eg— dressed ones (gold 
lines) used as input to obtain the self-energy. The doping 
dependence of the two coherent magnetic bands is in ex- 
cellent agreement with experiments [20] and earlier mean 
field calculations [lj| , and captures the incoherent weight 
(the UHB/LHBs) at higher energies seen experimentally. 
Note that the incoherent weight is concentrated near the 
top and bottom of the bare LDA bands, leading to a 
nearly doping independent UHB-LHB splitting. 

Further insight is provided by Fig. 3(a), which com- 
pares the fuU QP-GW DOS (blue hne) and the Sq- 
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FIG. 2: (color online) Spectral intensity as a function of cj 
along the high-symmetry lines for several dopings (at tem- 
perature r = 0). Blue to red color map gives the minimum 
to maximum intensity. In each panel, the gold lines represent 
the renormalized magnetic bands (Eo-dressed). 



dressed quasiparticle DOS (red line), normalized to the 
same peak height, for a representative doping x = 0.04. 
The good agreement between various computations over 
most of the energy range indicates the high degree of self- 
consistency in the self-energy. The dressed DOS shows 
four well-separated peaks. A clear leading edge gap of 
~0.3 eV (~1 eV at X = 0) can be seen. A gap persists 
at higher doping up to a; = 0.18 (Fig. 1), although it is 
obscured in the DOS by band overlap. The importance 
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FIG. 3; (color online) (a) The QP-GW DOS (blue lines) 
is compared with Eo-dressed DOS (red lines) calculated at 
T — 0. The green lines show the DOS without the ver- 
tex correction, (b) Solid lines give the imaginary part of 
the total self-energy, while the dashed lines give the corre- 
sponding charge contributions, (c) Doping dependence of self- 
consistent values of U and ZU is compared with earlier mean 
field results 14, 21]. (d) Renormalization factor Z decreases 
linearly with doping. 



of the vertex correction is illustrated by the green line 
in Fig. 3(a), which shows that setting F = 1 reduces the 
weight in the U/LHBs. 

Figure 3(b) shows how the calculated imaginary self- 
energy E" evolves with doping. The solid lines give the 
total S", with the corresponding dashed lines giving the 
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shown]. In the underdoped region the extra splitting at 
high energies in the self energy is related to spin-charge 
separation. The spin response is significant at low energy 
for all doping but the charge contribution is nearly zero in 
the low energy region for the lower doping and becomes 
finite only at higher energies above ~ 3 eV. As doping 
increases, the charge response moves toward the Fermi 
level and increases in contribution to its total value. At 
X = 18%, when the AFM gap vanishes, the charge and 
spin susceptibility become equal. Note that the broaden- 
ing of the self-energy evident in Fig. 3(b) is reflected in 
the increasing broadening of the Mott gap feature with 
doping. The shift of the peak in the imaginary part of 
the self-energy towards Ep in Fig. 3(b) reflects the dop- 
ing dependence of the MIR feature in Fig. 1. 

Fig. 3(c) describes the doping dependence of U. Al- 
though the renormalized Hubbard parameter ZU follows 
almost the same doping dependence as in mean field 
calculations [ij, |21|, the bare U displays considerably 
weaker doping dependence away from half-filling. The 
renormalization factor Z in Fig. 3(d) actually increases 
with underdoping. |22l| This can be readily understood. 
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FIG. 4: (color online) The computed DOS at various dopings 
are compared with the corresponding QMC results (blue lines) 
for X = 0.023] and for x = 0.05 to a; = 0.20 [H]. The red 
lines in each panel give our result for t' = 0, whereas the green 
lines in (b) and (c) are for t' = — 0.3t. 



our approximate self-energy calculation, the weight of the 
UHB is generally underestimated. 

In summary, we find that spin-wave dressing of the 
quasiparticles explains the incoherent U/LHB features 
seen in various experiments including the recent experi- 
ments observing waterfall effects in the cuprate spectra. 
The self-energy corrections not only renormalize the large 
widths of the LDA dispersions but also restore the resid- 
ual incoherent spectral weight associated with U/LHBs. 
In the underdoped regime, we show that the coherent in- 
gap bands reproduce both the four-band behavior seen in 
quantum cluster calculations and the magnetic gap col- 
lapse found in the mean-field calculations and a variety 
of experiments. The puzzling persistence of the Mott gap 
in optical spectra, even as the magnetic gap collapses, is 
thus reconciled. The fact that our calculations work so 
well confirms that the cuprates can be understood within 
the intermediate coupling regime, with U much less than 
twice the bandwidth. |25l| . 
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As we lower doping, the AFM gap increases leading to a 
decrease of the spectral weight near the Fermi level. This 
causes a reduction of the real part of the self-energy, shift- 
ing the peak in S' towards higher energy. The resulting 
slope decrease leads to a larger Z. 

Existing QMC calculations on the cuprates generally 
employ simpler hopping parameter sets, so we have re- 
peated our calculations with the same band parameters 
for a quantitative comparison in Fig. 4. All the QMC re- 
sults (blue lines) are obtained for t' — Q, U = St with 
a momentum dependent self-energy correction [23, |24| . 
We have obtained the corresponding DOSs (red lines) 
for i' = for S calculated at a fixed momentum of 
k = (7r/2,7r/2), using a doping dependent renormalized 
U (ZU) very similar to the one found for the electron 
doped case in Fig. 3(c). [Note that the QMC automati- 
cally generates a renormalized U] . Our result reproduces 
the QMC very well for cc = in Fig. 4(a) where a promi- 
nent electron-hole symmetry is observed. At higher dop- 
ing, the QMC results show a relatively smaller coherent 
peak above the Fermi level, whereas our result continues 
to exhibit near electron-hole symmetry in the in-gap re- 
gion as might be expected for i' = dispersion. We can 
mimic a momentum dependent self-energy by including 
t' — — 0.3i, and find that this provides better agreement 
with the QMC results as shown by the green lines in 
Figs. 4(b) and 4(c). Finally, at a; = 0.20, the pseudo- 
gap collapses and our result with t' = agrees very well 
with the QMC in Fig. 4(d). As might be expected from 
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